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Abstract. The aim of this article is to draw attention to 
the importance of the electric current loss in the energy 
output of radio pulsars. We remind that even the losses 
attributed to the magneto-dipole radiation of a pulsar in 
vacuum can be written as a result of Ampere force action 
of the electric currents flowing over the neutron star sur- 
face (Michel 1991, Beskin, Gurevich & Istomin 1993). It 
is this force that is responsible for the transfer of angu- 
lar momentum of a neutron star to an outgoing magneto- 
dipole wave. If a pulsar is surrounded by plasma, and there 
is no longitudinal current in its magnetosphere, there is 
no energy loss (Beskin, Gurevich & Istomin 1993, Mes- 
tel, Panagi & Shibata 1999). It is the longitudinal current 
closing within the pulsar polar cap that exerts the retar- 
dation torque acting on the neutron star. This torque can 
be determined if the structure of longitudinal current is 
known. Here we remind of the solution by Beskin, Gure- 
vich & Istomin (1993) and discuss the validity of such an 
assumption. Finally, it is shown that the behaviour of the 
recently observed "part-time job" pulsar B1931+24 can 
be naturally explained within the model of current loss 
while the magneto-dipole model faces difficulties. 



1. Introduction 

The recent observations of "part-time job" pulsars 
(Kramer et al 2006) such as pulsar B1931+24 has drawn 
attention to the particular mechanism of the energy losses 
of pulsars. In this article we summarize some results ob- 
tained for the model of current losses (Beskin, Gurevich 
& Istomin 1993) and the consequences of it. 

Radio pulsars are definitely nonaxisymmetric objects. 
However, the most results both in electrodynamics of the 
pulsar magnetosphere and in neutron star statistics were 
obtained under the assumption that the magnetic axis is 
parallel to the rotational one. Taking into account the 
inclination angle one can change qualitatively the conse- 
quences of the standard model. E.g., one can show that 
for the orthogonal rotator with the local GJ longitudinal 
electric current the light surface (where |E| 2 = |B| 2 ) must 
be located in the very vicinity of the light cylinder. In this 



case it is impossible to prolong the MHD flow up to in- 
finity, so the effective energy conversion and the current 
closure is to take place near the light surface. 

2. Magneto-dipole loss 

At first, the energy loss mechanism of radio pulsars has 
been connected with the magneto-dipole radiation (Pacini 
1967). Indeed, the magneto-dipole radiation power 



Wn 



1 BZ&R 6 2 

K 3 sm X, 



(1) 



where x is the angle between rotational and magnetic axis, 
R ~ 10 km is a neutron star radius, and is a pulsar 
angular velocity explains pulsar activity and observed en- 
ergy loss for expected large magnetic field near the surface 
B Q ~ 10 12 G. 

Let us recall that the physical reason of such en- 
ergy loss is the action of the torque exerted on the pul- 
sar by the Ampere force of the electric currents flowing 
over the neutron star surface (Istomin 2005). The elec- 
tric and magnetic fields in the outgoing magneto-dipole 
wave in vacuum can be found by solving the wave equa- 
tions V 2 B + n 2 /c 2 B = and V 2 E + fi 2 /c 2 E = with 
the boundary conditions stated as the fields corresponding 
components E t and B n being continuous through the neu- 
tron star surface. Inside the star one can consider magnetic 
field as homogeneous, and find the corresponding electric 
field using the frozen-in condition. As a result, the discon- 
tinuity of {B t } and {E n } give us the surface charge <r s 
and the surface current 



Js = ^[n,{B}]. 



The Ampere force exerts the torque 

K = - c Ar,[J 8 ,B]]dS 



(2) 



(3) 



on the neutron star. The energy loss of a pulsar due to 
this torque is equal to ((T|). Thus, it is the surface current 
that is responsible for the angular momentum transform 
from a neutron star to an outgoing magneto-dipole wave 
(Michel 1991, Beskin, Gurevich & Istomin 1993). 
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Thus, a pulsar in vacuum loses its rotational energy 
due to angular momentum transform to the electromag- 
netic wave at the rate given by (Q]). However, this is not so 
if the pulsar magnetosphere is filled with plasma and there 
is no longitudinal current in the magnetosphere. As was 
shown (Beskin, Gurevich & Istomin 1993, Mestel, Panagi 
& Shibata 1999), in this case the Poynting flux through 
the light cylinder is equal to zero. Indeed, as the ideal 
conductivity condition is applicable not only inside the 
neutron star but outside as well there is no magnetic field 
discontinuity at the star surface. Consequently, there is no 
Ampere force acting on a pulsar and, hence, there is no en- 
ergy loss. For zero longitudinal current the light cylinder 
is a natural boundary of the pulsar magnetosphere. 



3. Current loss 

In this section wc remind the exact solution for the surface 
current within the polar cap presented in the monograph 
by Beskin, Gurevich & Istomin (1993). As it was shown in 
the previous section, the neutron star retardation is due 
to Ampere force Fa = J s x B/c. If the magnetosphere is 
filled with plasma, the surface current J s is flowing within 
magnetic polar cap only. This surface current closes the 
volume longitudinal current in the magnetosphere and the 
return current flowing along the separatrix between open 
and closed field lines region. 

In order to write the equation for the surface current, 
the several assumptions must be made. We assume that 
the conductivity of the pulsar surface is uniform, and the 
electric field E s has a potential, so that the surface cur- 
rent can be written as J s = V£'. Using the stationary 
continuity equation divJ = where dJ z /dz is equal to 
the volume current i\\Bo flowing along the open field lines 
one can obtain 



The solution of homogeneous equation is 



V? 2) £' = -i\\Bo. 



Making the substitution x = sin 9 m and introducing 
the non-dimensional potential £ = 47r£' / BqR 2 H and cur- 
rent io = — 47rij|/ni? 2 we get 
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Here 6 m and ip m are polar and azimuth angles with respect 
to magnetic axis. 

Equation j5]) needs a boundary condition. This bound- 
ary condition results from the proposition that there is 
no surface current outside the magnetic polar cap. This 
means that 

£ [xo((p m ), ip m ] = const, (6) 

where xo(ip m ) is the polar cap boundary. Indeed, let us 
suppose that there is no such boundary condition for 
the potential £. In this case we get V? 2 )£ = r.h.s with 
r.h.s. = io inside the polar cap and r.h.s = outside it. 
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where / n (^) m ) = (a n cos nip m +b n sinmpm) with a n and b n 
being arbitrary constants. In this case the surface current 
J s oc V£ is circulating over the whole neutron star surface 
resulting in arbitrary energy loss. However, this means 
that there is a potential drop between different points of a 
neutron star surface which inevitably leads to the volume 
current in the region of closed field lines in the magneto- 
sphere. This contradicts to the assumption that there are 
no longitudinal currents flowing in the region of closed 
field lines. Thus, the boundary condition © is to be pos- 
tulated. In this case the jump in the potential derivative 
{V£} at x = xo(ipm) gives us the current flowing along 
the separatrix. As we see, it is defined uniquely by the 
longitudinal current in the region of open field lines and 
by condition that no longitudinal current can flow in the 
region of the closed field lines. 

For arbitrary inclination angle x the electric current 
io can be written as a sum of its symmetric is and anti- 
symmetric ix components. The anti-symmetric current 
begins playing the main role when the pulsar polar cap 
crosses the surface where the Goldreich- Julian charge den- 
sity pgj = — O • B/27rc changes sign. This condition can 
be written as 

7T [OR 

X=77-\/— ■ (9) 



For example, taking the Goldreich- Julian current density 



(4) lGj(x, ipm) 
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(and assuming xo = const) we obtain the following solu- 
tions of the Dirichlet problem (JH])-© for the symmetric 
and anti-symmetric volume currents respectively: 
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The torque exerted by the surface current over the 
neutron star can be written as 



(13) 



where Bo is the dipolc field near the neutron star surface. 
Let us expand the braking torque K over the orthogonal 
system of unit vectors e m , rii, and ri2. Here e m is a unit 
vector along the magnetic moment m; vector ri! is per- 
pendicular to the magnetic moment and lies in the plane 
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of the magnetic moment and the rotational axis; vector 
112 complements these to the right-hand triple: 



K = ^ + K±n 1 + A' t n 2 . 



(14) 



A-f plays no role in Euler equations that describe the ro- 
tational dynamics of the decelerating neutron star. As a 
result we have (Beskin, Gurevich & Istomin 1993) 
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and A = xcos ip m d^jdx — smtp m dt;/ dtp m . 

As one can find (Beskin, Gurevich & Istomin 1993), 
the torque component Ai is equal to zero equivalently 
for arbitrary shape of the polar cap due to the boundary 
condition ([6]). Thus, the values A|| and K± can be written 
as 
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Here the coefficients /i| and /jj^ depending on the shape of 
the polar cap are much less than unity, and the coefficients 
Cj| and cj_ ~ 1. 

We can now find the derivatives of the angular velocity 
Q and of the inclination angle x of a neutron star through 
the Euler dynamics equations: 



dn 

J r —r- = Am cos x + Ax sin y, 

dt 1 

dx 

J r fi— — = Aj^ cosx — A|| sinx- 



(21) 
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For the inclination angles x not too close to 90° (i.e., for 
cosx > (OA/c) 1 / 2 ) when the anti-symmetric current plays 
no role in the neutron star dynamics we obtain 
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As a result, for homogeneous current density within 
open magnetic field line region ig = j\\/jGJ — const where 
3gj = cpgj we have 



W c = -i- *s cosx- 



(25) 



Here /* is the non-dimensional area of a pulsar polar cap: 
•Scap = /*7r(fiA/c)A 2 . It depends on the structure of the 
magnetic field near the light cylinder. For a pure dipole 
magnetic field (and aligned rotator) /* = 1, and for a mag- 
netosphere containing no longitudinal currents /* changes 
from 1.592 for the aligned rotator, x = 0° (Michel 1991), 
to 1.96 for an orthogonal rotator, x = 90° (Beskin, Gure- 
vich & Istomin 1993). Recent numerical calculations for 
an axisymmetric magnctosphere with non-zero longitudi- 
nal electric current give /» « 1.23 - 1.27 (Gruzinov 2005, 
Komissarov 2006, Timokhin 2006). If the singular point 
separating open and close field lines can be located in- 
side the light cylinder, the value /* can be even ^ 1. As 
the Goldreich- Julian charge density near the polar cap is 
proportional to cosx, one can write 
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On the other hand, for x ~ 90° when the anti- 
symmetric current plays the leading role we obtain 
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As we see, the energy loss of the orthogonal rotator is 
VlR/c times smaller than that of the aligned rotator. 

4. PSR B1931+24 

The recent discovery of the "part-time job" pulsar PSR 
B1931+24 (Kramer et al 2006) with tl on /tl oS w 1.5 shows 
that the current loss is indeed playing an important role 
in the pulsar energy loss. If we assume that in the on- 
state the energy loss is connected with the current loss 
only and in the off-state with the magneto-dipole radiation 
(in which case the magnctosphere must be not filled with 
plasma) we get 
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It gives x ~ 60°. On the other hand, if we assume 
the Spitkovsky's relation for the on-state energy loss 
(Spitkovsky 2006) 



4 c 3 1 
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we obtain 



Oc 

a 



off 



3(1 + sin 2 x) 
2 



sin 2 x 



(29) 



(30) 



Clear, this ratio cannot be equal to 1.5 for any inclina- 
tion angle. This discrepancy can be connected with that 
fact that all the numerical calculations produced recently 
contain no restriction on the longitudinal electric cur- 
rent. As a result, current density can be much larger than 
Goldreich-Julian one. 
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5. On the magnitude of a surface current 

As we have shown, the current loss plays the major role 
in the pulsar dynamics. In particular, the behaviour of 
the pulsar B1931+24 can be naturally explained within 
this model. The current loss model have some important 
consequences: 

1. the energy loss of an orthogonal rotator is VlR/c times 
smaller than for an aligned rotator. This is connected 
with the boundary condition j6j that leads to almost 
full screening of the toroidal magnetic field in the open 
field lines region (see Beskin & Nokhrina 2004); 

2. consequently, during its evolution a pulsar inclination 
angle tends to 7r/2 where energy loss is minimal. 

On the other hand, it is known for the Michel's monopolc 
solution that in order to have the MHD flow up to infinity 
the toroidal magnetic field must be of the same order as 
the poloidal electric field on the light cylinder. If the lon- 
gitudinal current ju does not exceed by (ilR/c)^ 1 / 2 times 
the respective Goldreich- Julian current density (for the 
typical pulsars this factor approach the value of 10 2 ), the 
light surface |E| = |B| for the orthogonal rotator must be 
located in the vicinity of the light cylinder. In this case 
the effective energy conversion and the current closure are 
to take place in the boundary layer near the light surface 
(Beskin, Gurevich & Istomin 1993, Chiueh, Li & Begcl- 
man 1998, Beskin & Rafikov 2000). 

In order for these results being not true (for exam- 
ple, in order for the light surface being removed to in- 
finity) there must be a sufficient change in the current 
density value in the inner gap. We should emphasize that 
for the model with free particle escape it is hard to sup- 
port the current different than the Goldriech- Julian one. 
Indeed, since /?gj is the particle density needed to screen 
the longitudinal electric field, the current ju must be close 
jgj = epej- In order to change this value significantly 
we must support the plasma inflow in the inner gap re- 
gion (Lyubarskii 1992). For example, these particles can 
be produced in the outer gap. But it is obvious that for any 
poloidal field configuration the major number of field lines 
intersect the outer gap region inside the Alfvenic surface. 
Inside the Alfvenic surface the flow remains still highly 
magnetized. Thus, the deviation of current lines from the 
field lines is negligible. On the other hand, magnetized 
plasma can intersect the Alfvenic surface outwards only 
(sec Beskin 2006 for more detail). Thus, the outer gap can 
not significantly affect the current in the vicinity of the 
polar cap. 

Finally, it is necessary to stress that the recent numer- 
ical calculations (Gruzinov 2005, Komissarov 2006, Timo- 
khin 2006, Spitkovsky 2006) do not include into consider- 
ation that the longitudinal current density must be close 
to j'gj- In ah these works the authors assume that the cur- 
rent flowing through the cascade region can be arbitrary. 
However, if this is not so, and the current is indeed close 



to the Goldreich- Julian current, the structure of a magne- 
tosphere may be different from the ones obtained in the 
numerical simulations. 

6. Conclusion 

As we have seen, the current loss connected with the lon- 
gitudinal current flowing in the magnetospherc plays the 
main role in pulsar dynamics, and recent observations of 
"part-time job" pulsar supports this point. This evolution 
includes not only a neutron star retardation but also the 
sufficient change in the angle between magnetic and spin 
axis. We have seen as well that the model of current loss 
depends crucially on the distribution of the electric cur- 
rent and its value in the inner gap. For current loss model 
the inclination angle grows with time so a pulsar tends to 
be an orthogonal rotator. In this case the energy loss is to 
be VlR/c times smaller than for the aligned rotator. As a 
consequence, the light surface must be located in the very 
vicinity of the light cylinder. On the other hand, to realize 
the homogeneous MHD outflow up to infinity for the or- 
thogonal rotator the current density in the open field line 
region is to be much larger than j'gj ■ 
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